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ABSTRACT 

Based on SDSS DR8, we have selected a sample of 1728 red horizontal branch stars with |Z| < 3 kpc 
by using a color-metallicity relation and stellar parameters. The sample stars clearly trace a typical 
thick disk population with peaks at \Z\ = 1.26 kpc and [Fe/H] = —0.54. The vertical metallicity 
gradient of the thick disk is estimated in two ways. One is a fit to the Gaussian peaks of the 
metallicity histograms of the thick disk by subtracting minor contributions from the thin disk and 
the inner halo based on the Besangon Galaxy model. The resulting gradient is —0.12 ± 0.01 dex 
kpc^^ for 0.5 < |Z| < 3 kpc. The other method is to linearly fit the data based on stars with 
1 < 1^1 < 3 kpc being the main component of the thick disk. Five subgroups are then selected in 
different directions in the X — \Z\ plane to investigate the difference in the vertical metallicity gradient 
between the Galactocenter and anti-Galactocenter directions. We found that a vertical gradient of 
—0.22 ± 0.07 dex kpc^^ is detected for five directions except for one involving the pollution of stars 
from the bulge. The results indicate that the vertical gradient is dominant, but a radial gradient has 
a minor contribution for the thick disk population represented by RHB stars with 1 < |Z| < 3 kpc. 
The present work strongly suggests the existence of a metallicity gradient in the thick disk, which is 
thought to be negligible in most previous works in the literature. 

Subject headings: Stars:horizontal-branch - Starsdate-type - Galaxy: disk - Galaxy: kinematics and 
dynamics - Galaxy: structure 



1. INTRODUCTION 

The presence of the thick disk component in 
the Galaxy has been established for several decades 
(jGilmore fc ReidI Il983[ ). The age, metallicity and ro- 
tational velocity of thick disk stars have been investi- 
gated for various types of stars in the solar neighborhood. 
However, the formation mechanism of the thick disk is 
not well known. Actually, the metallicity gradient for 
the thick disk, being a basic input for the model of the 
Galactic disk, is not well studied in the literature. To 
the contrary, the traditional radial metallicity gradient 
of the Galactic thin disk is well investigated based on 
H// regions, Cepheids, OB stars, open clusters (young 
population) and planetary nebulae (old population). It 
is found that the radial gradients vary with time in the 
sense of a steeper slope in the early stages compared to 
the present value, and there is a change in slope at large 
Galactocentric distances i? > 10 kpc (jMaciel fc Costal 
120091) . 

Generally, the thick disk is typical for vertical distances 
from the Galactic plane of 0.5 to 1.5 kpc based on ob- 
servational data in the solar neighborhood. With a large 
spectroscopic survey of stars far away from the Sun, it 
is suggested t hat stars with 1 < jZI < 3 kpc dominate 
the thick disk (jAllende Prieto et al.ll2006l ). and the edge 
of the thick disk could reach up to \Z\ ~ 5.5 kpc. The 
large range in the \Z\ distribution enables us to inves- 
tigate the vertical metallicity gradient of the thick disk, 
which is not possible for the thin disk where most stars 
are located within \Z\ < 0.5 kpc. Such a study has re- 
cently begun and the existence of a metallicity gradi- 
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ent in the thick disk is a d iscrepant issue. For example, 
lAllende Prieto et"aLl ()2006[ ) proposed no metallicity gra- 
dient in the thick disk for 1 < |Z| < 3 kpc w ith a upper 
limit o n the gradient of 0.03 dex kpc~^, while iKatz et al.l 
()2011| ) investigated the metallicity distribution function 
(hereafter MDF) of the thin-thick-disk-halo system in 
two Galactic directions for several intervals between 
and 5 kpc above the Galactic plane and detected a ver- 
tical metallicity g radient of —0.068 ± 0.009 dex kpc~^. 
iSiegel et al.l (|2009[ ) suggested that the thick disk presents 
no metallicity gradient. But some models with a gradient 
are compatible with observations of the thick disk dom- 
inated b};_a^_20£ulation with 1 < \Z\ < 4 kpc as pointed 
out bv lKatz et al.l ()2011[ ). Meanwhile, in an indirect way, 
the mean metallicity for the thick disk varies with differ- 
ent heights from -0 48 ± 0.05 at 0.4 < \Z\ < 0.8 kpc 
(|Soubiran et al.|[2003l). to -0 685 ± 0.004 and -0.780 at 
1 < |Z| < 3 kpc (ISiegel et al.ll2009t ). which may indicate 
the existence of a metallicity gradient. Note that these 
works are mainly based on dwarf stars an d they have lim- 
itation s in some aspects. For example, iSoubiran et al.l 
(|2003D have high resolution but low signal-to-noise data 
for low heights of 0.4 < \Z\ < 0.8 kpc where both the 
thin and thick disks contrib ute s ignificantly to the MDF. 
lAllende Prieto et al.l ()2006l ) and llvezic eFall ()2008l ) were 
based on an SDSS DR6 sample where metallicity on the 
metal rich si de could be under es timated by .2-0.3 dex 
according to iBond et all (|2010f ). ISiegel et al.l (|2009[) are 
based on photometrically derived metallicities, which are 
not as good as spectros copically derived values. Re- 
cently, IKatz et alT (|2011| ) provide high quality data for 
the metallicity estimation, but the number of sample 
stars is limited to several hundreds. 
The aim of this work is to probe the metallicity gra- 
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dient of the thick disk population using red horizontal 
branch (hereafter RHB) stars selected from SDSS DR8. 
This study has some advantages over previous studies. 
First, the number of stars in the sample is on the order 
of thousands, and the metallicity is based on the up- 
dated version in SDSS DR8. Second, RHB stars are an 
interesting type of star, which are different from dwarf 
stars, with stellar distances easily obtained from their in- 
trinsic absolute magnitude. Their high brightness make 
them traceable to distances of up to 10 kpc. Finally, 
we investigate the metallicity distributions for thick disk 
stars in five different directions in the X — \Z\ plane. 
In particular, the comparison of the thick disk proper- 
ties between the Galactocenter and anti-Galactocenter 
directions has great potential for providing more strin- 
gent data for model comparisons. Thanks to the large 
survey area of the SDSS project, it is the first time that 
such properties in different directions can be studied and 
they will provide important clues for our understanding 
of the thick disk formation. 

2. THE SAMPLE STARS 

The selec t ion p rocedure of RHB stars is similar to 
Chen et all (|2010D . but uses S DSS DR8 (lAihara et alj 
201J rather than SDSS DR7 ijAbazaiian et al.ll2009n . 
Specifically, the first selection criterion is the color- 
metaUicity relatio n of (t? - r)n = 0.3 43(±0.039)[Fe/H] ^ 
0.829 presented in I Chen et al.l ()2010f ) with a deviation in 
{g — r)o to the above relation (5 < 0.20 mag. Then stars 
with (7o < 20 mag, [Fe/H] > —2.0 and signal-to-noise ra- 
tios of their spectra larger than 10 are used to obtain 
good quality data. Thirdly, RHB stars should have stel- 
lar parameters within the range of 4500-5900 K in tem- 
perature and 1.8-3.5 dex in logg. In the selection, pho- 
tometric data and stellar parameters (including metal- 
licities) in photometric and spectroscopic STAR tables 
of DR8 are adopted. Stellar distances of stars are cal- 
culated from g magnitudes and the absolute magnitude- 
me taUicity rela t ion o f Mg = 0.492[Fe/H] + 1.39 derived 
by IChen et al.l ()2009l ) with metallicities taken from the 
SDSS DR8 spectroscopic STAR table. The stellar dis- 
tance and position on the sky is then converted to the 
cartesian syst em of Galactic coord inates (X,Y,Z). 

As shown in lAihara et al.l ()2011t ). the photometric data 
and stellar parameters are updated in DR8 along with 
an enlarged survey area in the disk. Both improvements 
provide a better chance for us to probe the thick disk 
properties. With the aim to investigate the properties 
of the thick disk, the sample stars are selected to have 
|Z| < 6 kpc since the edge of the thick disk is at about 5.5 
kpc ab ove the Galactic plane as suggested bv lMaiewskil 
(119941 ). With this hmit, we extract a sample of 2,729 
RHB candidate stars from SDSS DR8. The left panel of 
Fig. 1 shows the distributions of distance to the Sun, ver- 
tical distance to the Galactic plane and Galactic latitude 
of the sample. In particular, the \Z\ distribution of RHB 
stars has a peak around \Z\ = 1.2 kpc and most of them 
are located within |5| = 6 — 16° with a weak tail extend- 
ing to \b\ = 40°. There is a hint of a minimum in the \Z\ 
distribution at \Z\ = 3, which is taken as the selection 
limit of the thick disk in many works. After that, \Z\ 
flattens out from \Z\ = 3 kpc to \Z\ = 6 kpc where we 
cut our sample. As we know that many previous works 
select thick disk stars with \Z\ < 3 kpc or \Z\ < 4 kpc, in 




Figure 1. Left: The distributions of distance, \Z\ magnitude and 
galactic latitude for RHB stars with \Z\ < 6 (the thick line) and 
\Z\ < 3 (the dashed line). Right: The go magnitude, {g — r)o color 
and proper motion distributions for stars with \Z\ < 3. 



the present work, we limit our sample to \Z\ < 3 kpc and 
1728 RHB stars are selected for studying the metallicity 
gradient. Most of them are limited to within 10 kpc and 
|6| < 20° as shown in the dashed lines in the left panel 
of Fig. 1. 

The right panel of Fig. 1 shows the go magnitude, 
(g — r)o color and proper motion (PM) distributions of 
our sample stars. In comparison with the selection cri- 
teria f or stars in the SDSS/SEGUE survev lYannv et al] 
(|2009f) . we found the main contribution of our sample 
stars comes from K (and partly red K) giants, which 
covers the metallicity range of [Fe/H] ^ —2.0 to the so- 
lar value (with l'^ > 0.07), has g magnitudes brighter 
than 18.5 mag, PM < 11 mas/yr and a color range 
of {g — r)o — 0.5 to 0.8 mag. These criteria will not 
produce a significant distance or metallicity bias in our 
sample selection, and thus we assume that these target 
selections in SDSS/SEGUE will not affect our result. 

3. ANALYSIS AND RESULTS 

3.1. The metallicity distribution of the thick disk 

The metallicity distribution of our sample stars is 
shown in the upper panel of Fig. 2. It shows that the con- 
tribution from the thin disk is significant with 550 stars 
having [Fe/H] > -0.25 and a peak at [Fc/H] = -0.12. 
The peak at \Z\ = 0.8 kpc in Fig. 1 partly corresponds 
to this population. The halo has a negligible contribu- 
tion with only 77 stars having [Fe/H] < —1.0. There are 
1101 stars with —1.0 < [Fe/H] < —0.25, corresponding 
to the thick disk, which have peaks at \Z\ = 1.26 kpc 
and [Fe/H] = -0.54. 
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Table 1 

Thick disk metallicity peak as a function of height 
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Figure 2. The [Fe/H] distribution of the sample stars before and 
after subtracting the contributions from the model thin disk and 
halo populations. 



In order to exclude the contribution of the thin disk 
and the halo, we obta in simulation data from the Be- 
sangon Galaxy model (|Robin et al.l 120031 ) using the se- 
lection criteria: (i) 12 < g < 20 mag; (ii) b = —30° to 
30° with a step size of 10°; (iii) / = 50° to 250° with a 
step size of 25°; (iv) absolute magnitude My = 0.0 — 1.2 
and spectral type from F to K; (v) population being ei- 
ther thin disk or halo; and (vi) the color-metallicity re- 
lation a nd stellar paramet er ranges for selecting RHB 
stars in iChen et al.l ()2010[ ). All these criteria fit our 
observed RHB star sample. Then, the metallicity dis- 
tributions of the simulated data for the thin disk and 
halo are normalized to the star numbers in the observed 
RHB sample for each population respectively. Coinci- 
dently, the simulated thin disk from the model has the 
same metallicity peak at [Fe/H] = —0.115 and the sim- 
ulated halo has negligible contribution to the metallicity 
distribution as shown by the dashed lines in the upper 
panel of Fig. 2. The MDF of the thick disk after sub- 
tracting the simulated thin disk and halo contributions 
is shown in the lower panel of Fig. 2. A Gaussian fit to 
the metallicity histogram is performed and has a peak at 
[Fe/H] — —0.54, which is quite close to the peak of the 
histogram. 

3.2. The metallicity gradient of the thick disk 

In our sample, there are significant numbers of stars 
at different \Z\ ranges of 0.5 — 1.0 kpc, 1.0 — 1.5 kpc, 
1.5 - 2.0 kpc, 2.0 - 2.5 kpc, and 2.5 - 3.0 kpc in the 
middle panel of Fig. 1. Thus, it provides a way to inves- 
tigate the metallicity gradient for < |Z| < 3 kpc. We 
perform the same procedure to each of the five differ- 
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ent \Z\ bins as that in Sect. 3.1 and obtain the Gaussian 
peak of the metallicity histogram for each \Z\ bin. Fig. 3 
shows the metallicity distribution for the 1.0 < \Z\ < 1.5 
kpc bin before and after subtracting contributions from 
the thin disk and halo components based on model data. 
The resulting peak metallicities from the Gaussian fits to 
the metallicity histograms for five different \Z\ bins are 
shown in Table 1, as well as the height range, the me- 
dian \Z\ and the star counts. The bin size for the plotted 
histograms are 0.05 dex and we have checked that a bin 
size of either 0.02 or 0.10 will not significantly change the 
Gaussian peak in the fit, but the distribution will deviate 
from a Gaussian function for one or two bins. The errors 
in the metallicity peak in Table 1 are estimated by the 
Gaussian dispersion of the metallicity peaks from 1000 
simulated samples randomly selected from the observed 
samples for each height bin using a bootstrap method. 
Fig. 4 shows the peak metallicity of the thick disk MDF 
as a function of median \Z\ and the metallicity gradi- 
ent is estimated from the weighted least squares fit to 
the data. The metallicity gradient is —0.12 ± 0.01 dex 
kpc~^ with an intercept of —0.34 dex. This gradient is 
slightly steeper than that of lKatz et al.l (j2011l ) who found 
-0.068±0.009 dex kpc~^ with an intercept of -0.46 dex. 

In order to estimate the errors from the Besangon 
Galaxy model, we vary the normalization factors in 
both the halo and the thin disk components by reduc- 
ing and enhancing the adopted value by 10%, which is a 
worst case scenario, and the modeled metallicity distri- 
butions deviate clearly from the observed distributions. 
Wc therefore obtain an additional four sets of Gaussian 
metallicity peaks of the thick disk for each bin. In total, 
we have five sets of metallicity distributions for five \Z\ 
bins, which produces 3125 metallicity gradients in the 
lower panel of Fig. 4. All these slopes lie within —0.11 
to —0.15, and thus the errors from uncertainties in the 
Besangon model do not have a significant effect on the 
metallicity gradients. 

3.3. The [Fe/H] versus \Z\ diagram in five different 
directions in the X — \Z\ plane 

IKatz et ahl ([2011) compared the metallicity distribu- 
tions in two different directions in order to investigate 
the origin of the gradient from either vertical or radial 
variations. They found that the vertical origin is strongly 
favored and the Kolmogorov-Smirnov test (hereafter KS 
test) shows similar metallicity distributions for two fields 
directed at M3 and M5 with diffe rent X — Z directi ons. 
Since both of the two directions in IKatz et al.l (j201lD are 
located in the anti-Galactocenter direction, it is inter- 
esting to make this comparison for different directions, 
including some Galactocenter directions. Meanwhile, the 
previous method to estimate the metallicity gradient is 
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Figure 3. The same as Fig. 2 but for the 1.0 < \Z\ < 1.5 bin. 
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Figure 4. The metallicity gradient for the thick disk for \ Z\ = 0—3 
kpc in this work (solid hne) as compared with that from iKatz et al.l 
|[20ni) (dashed line). 



based on the assumption that the thick disk does not 
include stars with high metaUicity ([Fe/H] ^ —0.12), 
which is arbitrarily assigned to be the thin disk popu- 
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Figure 5. Five star groups, Gl to G5, selected in the X — \Z\ 
plane. 



lation, and does not include stars with low metallicity 
([Fe/H] < —1.0), which is arbitrarily assigned to be the 
halo population. Actually, the thick disk can extend to 
a metallicity of [Fe/H] = 0.0 and a metal weak compo- 
nent of the thick disk is widely accepted to exist in the 
Galaxy. Therefore, we perform an alternative way to 
trace the metallicity gradient using the [Fe/H] versus \Z\ 
diagram, which is used in many works. In particular, it 
is important to know how the metallicity distributions 
vary between the Galactocenter and anti-Galactocenter 
directions since the scale length o f the thick disk typi- 
cally r anges from 2 .2 to 3.6 kpc ("M orrison et al.l [19901 : 
iRobin et al. 1996; Carollo et al. 2010). In our sample, it 
is possible to select stars in five different directions in or- 
der to investigate the properties of the thick disk on both 
sides of the solar circle. Fig. 5 shows the selection of five 
groups (Gl to G5) of stars in the X — \Z\ plane with the 
Gl and G2 fields toward the Galactic center, G4 and G5 
toward the anti-Galactocenter direction, and G3 is lim- 
ited to the solar region oi X — 8.3 to 9.3 kpc. Here X is 
along the line between the Sun and the Galactic center 
and is the absolute value of the vertical distance to 
the Galactic plane. 

Fig. 6 shows the [Fe/H] — \Z\ diagrams for all stars 
with 0.5 < \Z\ < 3.0 kpc in five directions. Stars with 
\Z\ < 0.5 kpc are excluded because they are missing in 
most directions and large contributions from the thin 
disk persist. The metallicity gradients are investigated 
by using mean metallicity at different height bins for five 
groups of stars. It is clear that all five groups show the 
existence of a metallicity gradient. A careful inspection 
of Fig. 6 shows the first two points with \Z\ < 1.0 kpc 
are quite different inside the solar circle (Gl and G2) 
compared to outside the solar circle (G4 and G5). In 
view of the fact that the thick disk is predominant for 
1 < ]Z] < 3 kpc while a star located between 0.5 and 1.0 
kpc has equal chances to be in either the thin or thick 
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Table 2 

The metallicity gradient in the formula of [Fc/H] = a \Z\ + b 
based on data for the five groups. 
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b 
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Figure 7. The five groups in the [Fe/H] — \Z\ plane. 
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Figure 6. The five subsamples in the [Fe/H] — \Z\ plane. 



disk populations according to the Besangon model of the 
Galaxy (Robin et al. 2003,), it is more reasonable to esti- 
mate the metallicity gradient among the different groups 
only for stars with 1 < |Z| < 3 kpc. Assuming that the 
variation of [Fe/H] with X can be neglected (see discus- 
sions in Sect. 3.4), the vertical metallicity gradients can 
be estimated with a linear fit to the data using the for- 
mula [Fe/H] = ajZj +5, which usually provides a higher 
slope than that from the peak of the distribution at each 
bin. For comparison, we derive such metallicity gradients 
for five groups, which is given in Table 2. The metallicity 
gradients for the four groups G2-G5 are similar with a 
slope of —0.22, but have decreasing intercepts from 0.02 
(G2) to —0.20 (G5). The regression to Gl gives a slope of 
—0.367 and an intercept of 0.07. A bootstrap method is 
used to estimate the errors of the slopes by adopting 1000 
simulated samples randomly selected from the observed 
five groups. The results are shown in Fig. 7. 

3.4. Vertical versus radial gradients 

When we fit linear relations to [Fe/H] as a function of 
\Z\, the variation of [Fe/H] with X is assumed to be neg- 



ligible. This may be the case for G2 to G5. Note that G2 
to G5 have different X variations: G3 has nearly a fixed 
X value and G2 has little X variation, while G4 and G5 
have significant X variation. Similar slopes in the [Fc/H] 
versus \Z\ fits for G2 to G5 indicate that the variation 
of [Fe/H] with X can be disregarded in the [Fe/H] ver- 
sus |Z| fits. This does not mean that there is no radial 
(hereafter X is referred to as radial) metallicity gradient. 
Instead, the decreasing intercept {b in Table 2) from G2 
to G5 indicates a negative radial gradient and a linear 
fit to the b versus median X for G2 to G5 gives a slope 
of —0.013 dex kpc~^. Such a small slope has no signifi- 
cant effect on the [Fe/H] versus \Z\ fits for our G2 to G5 
where X variations are less than 5 kpc. 

The small effect of X variation on the vertical metal- 
licity gradient for G2 to G5 can be further confirmed 
by similar metallicity distributions among these groups 
with different median X values for similar \Z\ ranges 
of 1 < \Z\ < 3 kpc. For this purpose, we investigate 
the cumulative distribution function (hereafter CDF) for 
[Fe/H] via the KS test among the five groups. For the two 
groups in the anti-Galactocenter direction G4 and G5, 
similar [Fe/H] CDFs wi th a jo- valu e of 0.1 4 are shown in 
Fig. 8. This agrees with lKatz et al.l (|201lD where the two 
directions toward M3 and M5 in the anti-Galactocenter 
region show similar [Fe/H] distributions. Since G4 is 
located inside of G5, similar [Fe/H] CDFs indicate the 
variation of metallicity with X is not detectable. Simi- 
larly, G2 and G3 have similar \Z\ distributions with little 
or no X variation within the group, and a KS test yields a 
p— value of 0.33 for the [Fe/H] CDFs based on stars with 
I < \Z\ < 3 kpc. Again, G2 is located inside of G3 and 
similar CDFs indicate X differences between G2 and G3 
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Figure 8. The cumulative distributions of G4 and G5 by KS test. Figure 9. Comparison of the [Fe/H] distribution between Gl and 

G5 by KS test for 1< \Z\ < 3 kpc. 



do not have a significant effect on their nietallicity distri- 
butions. Similar [Fe/H] CDFs are found for G3 and G4, 
but this is not the case for the comparison between Gl 
and G5. Fig. 9 shows the similar \Z\ but different [Fe/H] 
CDFs between Gl and G5, being the two extreme cases 
in opposite X directions. For —0.7 < [Fe/H] < —0.4, 
Gl and G5 have exactly the same CDFs while at both 
ends the differences in the [Fe/H] CDFs are significant. 
Two possible reasons may explain the special case of Gl. 
It has been suggested that the effect of [Fe/H] varia- 
tion with X cannot be neglected for a scale length of 
the thick disk between 2 and 4 kpc. Alternatively, it 
is expected that Gl has a high contribution of stars at 
both the metal poor end and the metal rich end from 
the bulge population with a wide metallicity range of 
— 2.0 < [Fe/H] < 0.5, which will steepen the gradient. In 
addition, it seems that Gl also has a special kinematic 
distribution in Fig. 10, which shows the radial velocity 
distributions for five groups with Gl having the largest 
radial velocity range and a peak significantly different 
from the other groups. Therefore, Gl is not included in 
this study. 

4. SUMMARY 

We have estimated the metallicity gradient for the 
thick disk population in two ways using RHB stars se- 
lected from SDSS DR8 data. The first method is based 
on the Gaussian peak of the metallicity distribution in 
five \Z\ bins by subtracting the contribution from the 
thin disk and halo via the Besangon Galaxy model. The 
slope is — 0.12±0.01 dex kpc~^ with an intercept of —0.34 
dex. The second method is linearly fitting the data in the 
[Fe/H] — \Z\ diagrams directly for stars from 1 < |^| < 3 
kpc, where the contribution from the thin disk and halo is 
not significant. Five groups in different directions can be 
separated to test for consistency, and we found that they 
give similar gradients of around —0.225 ±0.07 dex kpc~^ 
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Figure 10. Comparison of the radial velocity distribution for five 
groups at 1 < |Z| < 3 kpc. 



in four groups, with the only exception for the direction 
with —4 < X < 4. Both methods indicate the existence 
of a metallicity gradient in the thick disk. In our opinion, 
the first method gives a lower gradient because the model 
predictions suppress the existence of both metal rich and 
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metal poor components in the thick disk. We therefore 
favor the second method with a steeper gradient because 
it comes from the data directly and thus does not depend 
on the model assumptions. 
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